Impact of Point Source Clustering on Cosmological Parameters with CMB 

Anisotropics 
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The faint radio point sources that are unresolved in cosmic microwave background (CMB) 
anisotropy maps are likely to be a biased tracer of the large-scale structure dark matter distri- 
bution. While the shot-noise contribution to the angular power spectrum of unresolved radio point 
sources is included either when optimally constructing the CMB angular power spectrum, as with 
WMAP data, or when extracting cosmological parameters, we suggest that clustering part of the 
point source power spectrum should also be included. This is especially necessary at high frequen- 
cies above 150 GHz, where the clustering of far-IR sources is expected to dominate the shot-noise 
level of the angular power spectrum at tens of arcminute angular scales of both radio and sub-mm 
sources. We make an estimate of source clustering of unresolved radio sources in both WMAP and 
ACBAR, and marginalize over the amplitude of source clustering in each CMB data set when model 
fitting for cosmological parameters. For the combination of WMAP 5-year data and ACBAR, we 
find that the spectral index changes from the value of 0.963 ± 0.014 to 0.959 ± 0.014 (at 68% c.l.) 
when the clustering power spectrum of point sources is included in model fits. While we find that 
the differences are marginal with and without source clustering in current data, it may be necessary 
to account for source clustering with future datasets such as Planck, especially to properly model 
fit anisotropies at arcminute angular scales. If clustering is not accounted and point sources are 
modeled with a shot-noise only out to / ~ 2000, the spectral index will be biased by about 1.5cr. 

PACS numbers: 98.70.Vc,98.65.Dx,95.85.Sz,98.80.Cq,98.80.Es 



I. INTRODUCTION 

As discussed in a variety of papers, unresolved radio 
point sources are an important foreground in temper- 
ature anisotropy maps of the cosmic microwave back- 
ground (CMB) [H, 0, i, g. With Wilkinson Microwave 
Anisotropy Probe (WMAP) data [f|, the difference in 
the CMB power spectra determined at various frequency 
channels and the cross power spectra between the chan- 
nels, has allowed the unresolved point source contam- 
ination to be constrained with an amplitude ^4 ps = 
(0.011 ± 0.001)^K 2 -sr for the foreground power spec- 
trum, when scaled to the Q-band. In the WMAP analy- 
sis, this point source amplitude is taken to be a constant 
in Ci , similar to the case of a shot-noise type power spec- 
trum for unresolved point sources. This shot-noise, with 
an appropriate scaling in frequency, is then removed from 
each of the power spectra when constructing the final 
WMAP temperature anisotropy power spectrum Q- 

While the WMAP estimate on the point source cor- 
rection is consistent with a point source power spectrum 
dominated by the shot-noise, this estimate is dominated 
by measurements relative to the Q-band [||. Since the 
WMAP temperature anisotropy power spectrum is based 
on V- and W-band data Q , and the point source correc- 
tion has a larger uncertainty in V- and W- bands Q, it 
is unclear if a simple shot-noise correction fully describes 
point sources in the WMAP temperature anisotropy 
power spectrum. To account for uncertainty in the am- 
plitude of point-source shot-noise in parameter estimates, 
the WMAP likelihood contains an additional marginal- 



ization of the uncertainty of ^4 ps , but the best-fit ampli- 
tude of point-sources remain fixed to the a priori deter- 
mined value [3] • We also note that alternative approaches 
have been considered to estimate the amplitude of point- 
sources 0, [T3| , though these works also concentrated on 
establishing the shot-noise correction. 

Beyond the shot-noise, unresolved radio point sources 
are likely to have a clustered distribution on the sky as 
they are expected to be a biased tracer of the large-scale 
structure. Thus, the angular power spectrum of sources 
contains not just a shot-noise but also a clustering piece 
determined by the dark matter power spectrum, point 
source bias, and the redshift distribution. Existing calcu- 
lations suggest that the shot-noise part of the power spec- 
trum from bright, rare sources dominates clustering at 
low radio frequencies, especially when the flux threshold 
for point source removal is at the level of ~ 1 Jy @, |Tl[ • 
Thus, the assumption of a shot-noise point source con- 
tribution to the angular power spectrum of temperature 
anisotropies is likely to be adequate for low-frequency 
bands of WMAP such as the Q-band, but may not be 
appropriate at high frequencies, such as the W band, 
whose data are used in the temperature power spectrum. 

The approach using a shot-noise spectrum with an un- 
certainty that is marginalized over in the WMAP likeli- 
hood is is bit different from the approach advocated by 
the WMAP team to account for another foreground in 
CMB data involving the Sunyaev-Zel'dovich (SZ) effect 
from galaxy clusters. There, the angular power spectrum 
is estimated based on a model for the cluster distribution 
and gas properties [l2| , with an overall uncertainty in the 
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amplitude of the SZ power spectrum captured by a free 
parameter which is then freely varied as a nuisance pa- 
rameter when best-fit cosmological parameter values and 
their uncertainties using a Markov-Chain Monte-Carlo 
(MCMC) code [llllllj. 

Since the clustering component of the angular power 
spectrum of unresolved radio sources may be important, 
it could be that simply including point-sources as a shot- 
noise correction in the V/W-band WMAP temperature 
anisotropy power spectrum results in biased estimates of 
cosmological parameters, especially for parameters like 
the spectral index of density perturbations, which has 
been discussed previously in the context of uncertainties 
related to point-source shot-noise amplitude [9(. More- 
over, current CMB analyses make use of the combina- 
tion of datasets such as WMAP and ACBAR which have 
different treatments related to how point sources are ac- 
counted in the parameter fits. At the ACBAR frequency 
of 150 GHz [TH, the clustering of sources may need to be 
included properly, especially given that the high angular 
resolution of ACBAR also allows removal of sources down 
to a lower flux density level, where the shot-noise asso- 
ciated with rare, bright sources may be subdominant. 
The clustering of point sources could also account for 
some fraction of the excess arcminute-scale anisotropics 
detected by ACBAR [II [H. 

Given the lack of adequate details related to the exact 
clustering power spectrum of radio sources in datasets 
such as WMAP and ACBAR, we make a general estimate 
of the angular power spectrum of radio point-sources and 
include the overall amplitude of clustering as an addi- 
tional nuisance parameter to be included and marginal- 
ized over when estimating for cosmological parameters. 
For this, we make use of number counts at 95 GHz [l8[ 
and assume the redshift distribution of high-frequency 
radio sources follows the same distribution as estimated 
for NVSS sources at low frequencies [2(j. As some of 
these assumptions are likely to be invalid to some extent, 
we do not fix the clustering spectrum to our model but 
allow the overall amplitude to vary and marginalize over 
that uncertainty when constraining cosmological param- 
eter. Thus, the uncertainty in our predictions related to 
clustering of sources is unlikely to dominate and we con- 
firm this by noting that the differences to best-fit param- 
eter values with point source clustering included are not 
significant, especially for the case of combined WMAP 
and ACBAR data. Moreover, our estimate of cluster- 
ing is consistent with the allowed level of point source 
correction in V- and W-band combination of WMAP, as 
measured in terms of differences in the power spectra @ . 

While the differences in estimated cosmological param- 
eters are small, a proper estimate of cosmological param- 
eters with clustering included is useful for a proper sta- 
tistical analysis on important scientific results such as on 
the extent to which the spectral index of density pertur- 
bations n s is different from the Harrison-Zcl'dovich value 
at -1. While the impact on current datasets is small, for 
future data such as Planck that probe down to smaller 



scales over a wide range of frequencies, we suggest that it 
will be necessary to account for clustering of point sources 
when model fitting cosmological parameters. 

This paper is organized as follows: we first discuss 
the angular clustering power spectrum of radio points 
sources. Section III discusses model fits to recent CMB 
anisotropy data from WMAP and ACBAR. We discuss 
our results and conclude with a summary in Section IV. 



II. CLUSTERING OF RADIO POINT SOURCES 

The angular power spectrum of radio sources, in units 
of (/j,K) 2 , generally contains two components [2ll |22| 
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where wi is the Legendre transform of the angular corre- 
lation function w(9) of unresolved radio point sources, I 
is the average intensity (in flux units) produced by these 
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and the conversion factor from flux to antenna temper- 
ature using the CMB black-body spectrum, B„(T — 
2.726K), is 
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where x = hv/ksT = z^/56.84 GHz is the dimension- 
less frequency. This conversion can be simplified as 
dB v jdT = [{99.27Jysr- 1 )/nK]x 4 e x /{e x - l) 2 . 

Since we know little about the clustering of radio 
sources below the point source detection limit W- and 
V- bands of WMAP and at 150 GHz of ACBAR, wc 
make use of a simplified set of assumptions to estimate 
the source clustering. In general, the angular power spec- 
trum of the source sources can be written with the halo 
model 1231 such that wi is 



dr 2 n 2 {z) 
dZ Tz a {z) ^(z) Ps - 
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(4) 



where P ss (k, z) is the three-dimensional power spectrum 
of radio sources as a function of redshift. In the halo 
model, source clustering at large angular scales can be 
described with the linear matter power spectrum scaled 
by a constant and scale-free bias factor: 



P S8 {k) « blP lia (k) , 



(5) 



where the source bias factor, when combined with an es- 
timate of the number density of sources, provide some 
information on the halo mass scale associated with those 
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sources through the luminosity- or flux-averaged halo oc- 
cupation number (N(M,z)), halo bias b^ioiM, z), and 
the halo mass function dn/dM (23J: 

bs = J_ f dM-^-(z) 6 halo (M, z)(N(M, z)) . (6) 
n g J dM 

At small angular scales, clustering traces the non-linear 
power spectrum generated by the so-called 1-halo term. 
Separating the occupation number to central and satellite 
radio sources, {N{M)) = (N s ) + (N c ), the 1-halo power 
spectrum is 



P ih (k) 



dM n(M) 



(7) 



2(N s )(N c )u(k\M) + {N s ) 2 u 2 (k\M) 



Here, u(k\M) is the normalized density profile in Fourier 
space. At deeply non-linear scales, however, the shot- 
noise term is expected to dominate the clustering spec- 
trum, but the transition scale may lie larger than the 
shot-nose amplitude. The above form of the 1-halo term 
allows us to easily understand a simple behavior. If radio 
sources occupy dark matter halos such that there is only 
one source per halo, regardless of the halo mass, then 
with N s = 0, P lh = 0. Thus, the 1-halo term only exists 
to the extent that more than one radio source occupies a 
halo. While there is limited information on the halo oc- 
cupation properties of radio sources at the frequencies of 
interest, observations at frequencies around 30 GHz sug- 
gest that multiple radio sources are found in large dark 
matter halos such as groups and clusters, though at 30 
GHz, the central galaxy tends to be the dominant bright 
source in most galaxy clusters [24| . 

Given the lack of detailed knowledge on the clustering 
of radio sources or even ingredients such as luminosity 
functions or exact redshift distributions that can be used 
to generate a reliable halo model for the radio source 
population using approaches such as the conditional- 
luminosity functions that are used to describe clustering 
of optical or IR and far-IR galaxies [25[, we make sev- 
eral approximations. First we note that at large angu- 
lar scales, Ci « I 2 {b 2 )w\ in [21(. To calculate the angular 



power spectrum, we assume that unresolved sources trace 
the same large-scale structure as low-frequency NVSS 
sources and estimate / by integrating over the number 
counts at 95 GHz as estimated by [18|]- We use 95 GHz 
as a first estimate here since it is close to both WMAP 
channels on one end and ACBAR at the other end. We 
make use of the redshift distribution estimates for NVSS 
to calculate clustering at high frequencies [2(|. In ad- 
dition to linear clustering, we also include a non-linear 
correction to the angular clustering using a 1-halo model 
that assumes a simple power-law occupation number for 
satellite galaxies with N a (M) ~ M 13 with (3 = 0.85 when 
M > 10 12 - 5 M . The typical bias factor estimate for 
sources from this occupation number is about 1 at z ~ 1. 
Before calculating anisotropics for CMB, we verified that 



our prediction for source clustering, when applied for low- 
frequency sources, generally agrees with measurements 
from the literature [26|, [13] • 

In Fig. 1, we show the angular power spectrum of radio 
sources as fluctuations in the CMB temperature C; , and 
a comparison to the difference in power spectra of V and 
W-bands of WMAP @. For this comparison, we follow 
the same procedure as the WMAP analysis Q and scaled 
the power spectrum to Q-band (40.7 GHz) and estimate 
^ ps = HQy 2 CP(Q) with r{v) = (e x - l) 2 /(x 2 e x ) with 
a numerical value for the Q-band of r(Q) = 1.089. In 
Fig 1, we have converted our estimate of C;'s from 95 
GHz counts to Q-band with average spectral index of 
bright resolved WMAP sources with a f=a —0.09 [sj], with 
the scaling v a ~ 2 r(y) for temperature units instead of in- 
tensity. In addition to clustering we also include the V- 
and W-band combined estimate of shot-noise in WMAP 
data with a value of 0.007 /iK 2 -sr, when scaled to the 
Q-band. 

As shown in Fig. 1, the sum of this shot-noise and 
the clustering we estimate is consistent with the allowed 
amplitude of point source correction from Nolta et al. 
[(| . Though there are larger uncertainties in point source 
estimates of the V and W-band, we cannot simply rule 
out that the unresolved sources only contribute with a 
shot-noise type power spectrum. If Q-band is also in- 
cluded, as shown in Ref. [||, the differences are more 
consistent with a shot-noise power spectrum, as expected 
for low-frequencies since Q-band dominates such an es- 
timate. Since the final WMAP power spectrum is com- 
posed of V- and W-band data, we find that there is some 
motivation to include source clustering when estimating 
cosmological parameters. 

While the clustering amplitude is unconstrained, this is 
of little concern to us in the cosmological parameter esti- 
mation since we will not fix the clustering of radio sources 
to a pre-determined model, but parameterize the overall 
amplitude of clustering with a free parameter. Thus, 
when model fitting the data we parameterize the cluster- 
ing part as -Pwmap Cz and consider Pwmap as a nuisance 
parameter that captures all uncertainties in our calcula- 
tion, which includes the spectral index of sources from 
counts at 95 GHz to WMAP band, their redshift dis- 
tribution, parameters of the halo model, among others. 
When quoting cosmological parameter measurements, we 
marginalize the likelihood over Pwmap • This approach is 
consistent with how the WMAP team included the effect 
of SZ angular power spectrum in parameter estimation 
with a parameter ^4sz that is freely varied. Note that 
we only include a model for the source clustering since 
the CMB power spectrum released by the WMAP team 
already has a shot-noise removed from the data for point 
sources when combining V and W-data to a final power 
spectrum. 

In addition to clustering of sources as relevant for 
WMAP, we also include clustering of point sources as 
related to ACBAR data [l^|. In parameter estimation, 
unlike the WMAP team that fitted and removed a con- 
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FIG. 1: The amplitude of point source correction to the V- 
and W-band WMAP angular power spectrum. The data 
points show the measurement from the WMAP team scaled 
to the Q-band. We ignore the corrections with Q-band as the 
final power spectrum from the WMAP team uses only V- and 
W-band data. The dotted line shows our model for the point 
source clustering (scaled to Q-band) with Pwmap = 1, while 
the solid line shows the total clustering spectrum arising from 
point sources with clustering and the shot-noise. The shot- 
noise is taken to be the same as estimated by the WMAP team 
for V and W-band data with a value Aps = 0.007 /iK 2 -sr 



stant shot-noise spectrum for unresolved radio sources 
when estimating an optimal power spectrum from data, 
the ACB AR team included the shot-noise of radio sources 
as an extra component in their model fits. Thus, while 
we only include clustering spectrum of radio sources for 
WMAP, for ACBAR data, we include both a clustering 
spectrum and a shot-noise for radio sources. The shot- 
noise was taken to be consistent with estimates made 
by the ACBAR team and the clustering component was 
taken by simply frequency scaling the same WMAP spec- 
trum to 150 GHz with the same scaling as the one in- 
volved with the shot-noise part. 



III. COSMOLOGICAL PARAMETERS WITH 
CLUSTERED POINT SOURCES 

The method we use to estimate cosmological param- 
eters is based on the publicly available Markov Chain 
Monte Carlo package CosmoMC [l3| with a convergence 
diagnostics based on the Gelman and Rubin statistic. 
We used WMAP 5-year data [Hj (both temperature and 
temperature-polarization cross-correlation) alone and in 
combination with ACBAR data (la [. We only account 
for point sources in temperature anisotropics. Since 
WMAP polarization data do not probe small angular 
scales, where polarized point sources contribute, ignor- 
ing point sources in polarization is a safe assumption. 

In our estimates we make use of the flat ACDM cos- 



mological model with 6 cosmological parameters: baryon 
density f2f,/i 2 , dark matter density fl c h 2 , reionization op- 
tical depth r, ratio of the sound horizon to the angular 
diameter distance at the decoupling measured by 9, am- 
plitude of the curvature perturbation A s (with flat prior 
on log(A s )) and spectral index n s ; these two last param- 
eters are both defined at the pivot scale ko = 0.002/ Mpc 
as in ||. To this set we include Asz, the amplitude of SZ 
contribution, and two parameters Pwmap and Pacbar 
for the amplitude of point-source clustering. To study 
the impact of point sources on running of the spectral 
index and estimates of the tensor-to-scalar ratio, we also 
consider additional runs where these quantities are var- 
ied. 



A. WMAP and ACBAR data 

When estimating parameters with existing WMAP 
and ACBAR data, with point sources and SZ included, 
the total CMB anisotropy spectrum is 



Qtot _ £«CMB 



C, PS 



(8) 



The point-source angular power spectrum contains two 
parts as discussed: Cf s = Cf n + Cf, but since the 
WMAP team removed the shot-noise when combining 
data to a single estimate of the power spectrum, we take 



pWMAP _ (jCMB 



Cf + Cf z . There is a slight compli- 



cation here since Cf is a combination of the clustering 
in V- and W-bands, and we make the simple assump- 
tion here that the clustering of sources between these 
two bands can be scaled by a constant while the shape 
remains the same. The uncertainty in the variation of the 
point source clustering with frequency, to some extent, is 
not expected to be a significant issue since we allow the 
overall amplitude to vary with Pwmap- Note that the 
same complication exists for Cf z , but in this case the 
frequency dependence is known exactly. 

As described, in addition to WMAP 5-year data, we 
also include ACBAR data at large multipolcs. To avoid 
complicating the analysis when different datasets over- 
lap, which requires a calculation of the covariance ma- 
trix between different experiments, as they observe the 
same sky, in the likelihood calculation, we take the same 
approach as the WMAP team and use WMAP data out 
to I < 900 and ACBAR data from 900 < I < 2000. 
For ACBAR data, we make a separate estimate of the 
angular power spectrum by scaling the flux-cut of un- 
resolved point sources to be at the lower flux threshold 
and in agreement with previous shot- noise estimates [Pol ] . 
Again, we include an overall uncertainty in the ACBAR 
angular power spectrum of radio sources, in this case 
the sum of clustering and shot-noise terms of the power 
spectrum, with Pacbar- In Fig- 2, we show the angular 
power spectrum of CMB anisotropies with best-fit cos- 
mological model for WMAP and ACBAR data, as well 
as the two input power spectra for point source clustering 
with both Pwmap = Pacbar = 1. 
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TABLE I: Mean values and marginalized 68% c.l. limits for several cosmological parameters from WMAP and WMAP+ACBAR, 
with and without clustering of point sources (see text for details). 



Parameter 


WMAP 5-yr 


WMAP 5-yr 
with clustering 

Cf X PwMAP 


WMAP 5-yr 
with clustering 

Cf X (0 < PwMAP < 1) 


WMAP±ACBAR 


WMAP±ACBAR 
with clustering 
Cf x (0 < Pwmap < 1) 
Cf x (0 < Pacbar < 1) 


fl b h 2 


0.02277 ± 0.00062 


0.02397tHow4 


02S66+ 00084 

U.UZOUU_q .00083 


02269+ ao0059 
u.uiiuy_ .00060 


n n99Q«+ 000f3 3 

u.uiiyo_Q 00065 


Q c h 2 


0-1093+^63 


O.W25+° Zl 


0.1041 ±0.0069 


0.1103 ± 0.0059 


0.1092+°^ 


n A 


0.744ir 3 2 ° 9 


0.780 ± 0.034 


n 779 +0.031 
' ' ^—0.030 


0.740l°:°oi 


0.747l°;°f 7 


n s 


0.965 ±0.014 


0.949l°;°i 8 7 


0.953 ±0.016 


0.963 ±0.014 


0.959 ±0.014 


T 


0.087 ±0.017 


0.088 ± 0.017 


O-OSSlHir 


0.087l°;S 6 7 


0.086l 9 ;S 8 6 


A 2 


(2.39 ±0.10) • 10" 9 


(2.33 ±0.10) • 10" 9 


(2.35 ±0.10) ■ 10" 9 


(2.41 ±0.10) ■ 10~ 9 


(2.40 ±0.10) ■ 10~ 9 


o~s 


0.793 ± 0.036 


0.726 ± 0.056 


0.742 ± 0.047 


0.798 ± 0.033 


0.784 ± 0.033 




256+ 029 

U.ZOU_ 030 


0.22 ± 0.034 


n OOO+0-031 

U.ZZO_ .030 


2 60+ 0029 

U.ZDU_ .028 


n 9 ro+0.0027 
U.^OO_o.o28 


H 


72 1 + 2 - 7 
' z - 1 -2.6 


76.31" 


75.3 ± 3.4 


71.7 ±2.5 


72.5 ± 2.6 


Zreion 


11.0 ± 1.4 


10.5111 


10.7 ± 1.4 


11.0 ± 1.4 


10.8 ± 1.4 


to 


13.68±° ;S 


13.49 ±0.19 


13.54 ±0.16 


13.69I° ;« 


13.65 ±0.13 


A sz 


1-041^1 


1.00 ±0.68 


1.00 ± 0.68 


o.98lg;S 


n 01 +0-68 
U.yi_ .65 


-PwMAP 




< 1.38(2o-) 


< 0.93(2<j) 




< 0.46(2a) 


-Pacbar 










< 0.56(2a) 




10 1 00 1 000 



Multipole I 

FIG. 2: Best-fit CMB angular power spectrum for WMAP 5- 
year data and ACBAR with a comparison to measurements. 
We also show the input power spectra of point source cluster- 
ing for WMAP (middle line) and ACBAR data (bottom line) 
with Pwmap = Pacbar = 1 (see text for details) . In the case 
of WMAP, we show custering Cf part only as the shot-noise is 
removed from the data, while for ACBAR we show the total. 



Since wc only include the clustering term of unresolved 
point sources for WMAP data, we follow the WMAP 
team's approach on the shot-noise term and marginalize 
the likelihood over the uncertainty related to point source 
shot-noise term A ps using the public WMAP likelihood 
routine. This uncertainty only makes a small difference in 
best- fit cosmological parameters jf|. The results related 
to ACDM runs arc summarized in Table I. In the case 



TABLE II: Constraints on spectral index n s , running of 
the spectral index dn s /d\nk and tensor to scalar ratio r 
from WMAP±ACBAR with and without point source signal. 
These values are evaluated at k = 0.002 Mpc -1 . 



parameter WMAP±ACBAR WMAP± ACBAR 

with clustering 


n.s 

dris/dln k 


1.036 ±0.046 
-0.038 ±0.024 


1.041 ± 0.050 
-0.044 ± 0.025 


n s 
r 


0.981 ±0.019 
< 0.36(2a) 


+0.020 

u.y( '_o.oi9 
< 0.39(2cr) 



where we do not consider clustering of point sources, we 
essentially recover the same results as Ref . Q , with small 
differences at the level of O.ler, which we believe is due 
to differences in the numerical codes and the convergence 
criteria. 

With clustering of point sources included, however, the 
spectral index estimated with WMAP 5-ycar data alone 
changes from 0.965 ± 0.014 with point source shot-noise 
only to 0.949 ± 0.018 with source clustering in addition 
to the shot-noise from the WMAP likelihood. This is a 
difference of about la, but this large difference primarily 
comes from the fact that Pwmap is largely unconstrained 
by the data with a 2a upper limit of 1.38. The change 
in n s is captured by a similar change in as with values 
changing from 0.793 ±0.036 without clustering to 0.726± 
0.056 with clustering. If we put a prior that Pwmap is 
uniform between and 1, n s = 0.953 ± 0.016 and the 
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FIG. 3: Marginalized parameter constraints for WMAP without clustering (red line), WMAP with point source clustering 
(green line), and WMAP+ACBAR with point source clustering signal for both (blue line). From left to right, each of the 
panels show the constraint on the baryon density, cold dark matter density, the ratio of sound horizon to the angular diameter 
distance at the decoupling (top panels), optical depth, spectral index, amplitude of curvature perturbations (middle panels), 
and SZ normalization, WMAP point source normalization, and ACBAR point source normalization (lower panels), respectively. 
The spectral index and the amplitude of perturbations is measured at k — 0.002 Mpc -1 . 




FIG. 4: Two-dimensional marginalized distributions showing the 68% and 95% confidence level contours for n a vs. h 2 (left) 
and n a vs. as (right) with WMAP 5-years data alone (empty contours), WMAP 5-years data with clustered point sources (red 
contours), and WMAP 5-years data+ ACBAR data with point source clustering (blue contours). 
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difference from the case with clustering ignored is about 
- 0.8ct. 

With the addition of ACBAR data, the clustering am- 
plitudes of unresolved point sources in both WMAP and 
ACBAR are better constrained. Though we only use 
WMAP point source clustering model for WMAP data 
only and a separate model for ACBAR point sources 
at I > 900, both parameters are better constrained be- 
cause the combination of WMAP and ACBAR data pin 
down the overall cosmological model leaving less room for 
the point source piece to change in amplitude. In com- 
bination, with clustering of point sources included for 
both WMAP and ACBAR, we find n s = 0.959 ± 0.014, 
which is different from the WMAP+ACBAR value of 
n s = 0.963 ± 0.014 by about 0.3ct. 

In Fig. 3 we summarize likelihoods of the parameters 
involved. As shown there, when clustering is included 
the large differences on cosmological parameter estimates 
with WMAP data alone appear in n s ,Q,bh 2 and Q, c h 2 . 
However, as discussed for n s , once we include ACBAR 
data and with clustering of point sources both for WMAP 
and ACBAR, the probability distributions are more con- 
sistent with the WMAP data alone, but with clustering 
ignored. While it seems like large multipole data from an 
experiment such as ACBAR do not improve cosmological 
parameters estimated from WMAP, in our case, we do see 
an improvement by constraining the point source cluster- 
ing amplitude better. The associated contour plots for 
parameters that are mostly affected by clustering of point 
sources are summarized in Fig. 4 for combinations of n s 
vs. ilh/i 2 (left) and n s vs. erg (right). 

In addition to standard ACDM runs with a power-law 
power spectrum for density perturbations, we also study 
the impact of point sources on the running of the spectral 
index and on the tensor-to-scalar ratio, in addition to 
the main parameters of the ACDM cosmological model 
outlined in Table I. In Table II we summarize our results. 
Our results for the combination of WMAP and ACBAR 
without clustering are generally consistent with previous 
results [28| , but with minor differences such as a 2a upper 
limit on r of 0.36 instead of 0.4. The differences between 
with and without clustering are also minor and this is 
primarily due to the fact that we take the combination 
of WMAP and ACBAR. In Fig. 5 we summarize these 
results in contour plots with n s vs. running (dn s /dlnfc, 
left) and n s vs. r (right). 

While we find differences at the level of la for WMAP 
data alone with clustering of point sources, our results 
show that the differences are smaller and insignificant 
once WMAP data are combined with ACBAR data and 
using two separate estimates for point source clustering 
in WMAP and ACBAR data. In future, Planck data will 
observe CMB anisotropics down to smaller angular scales 
and extending to higher frequencies where clustering of 
sources becomes increasingly important [IJ HH ■ I n this 
case, it is clear that a simplified approach with a shot- 
noisc for unresolved point sources in Planck data may not 
be appropriate when extracting cosmological parameters. 



B. Planck mock data 

To understand how clustering of point sources impact 
cosmological parameter determination, we created sev- 
eral mock datasets with noise properties consistent with 
Planck 143 GHz channel of HFI and assuming the best- 
fit WMAP5 parameters [28| for cosmology. We model 
the point source clustering and the shot-noise by mak- 
ing use of existing high-frequency data as we did for 
ACBAR. While we only consider a single clustering spec- 
trum, at high-frequencies of Planck HFI, two separate 
populations of point sources are expected: radio, domi- 
nating at low- frequencies, and sub-mm or far-IR sources 
at high-frequencies [25| . Here, as we only have total num- 
ber counts at 150 GHz, without any information on hwo 
to separate the counts to the two populations, we make 
use of a single clustering spectrum. It will be necessary 
to return to this topic later once Planck data become 
available with additional information, from Planck and 
Herschel, on the far-IR population in HFI channels. 

We summarize our results related to Planck data in 
Fig. 6. In addition to the analytical model of point source 
clustering used in this paper based on the halo model, we 
also made use of publicly available Planck source maps 1 
from the Planck Working Sub-Group for Compact Source 
Fields to measure the angular power spectrum of points 
sources at the Planck HFI 143 GHz channel. These maps 
are derived from a model based on GallCS model 2 using 
the Mock Map Facility (MoMaF,[10|). The power spec- 
trum is computed after removing point sources with a 
flux greater than 72 mJy, 5a detection level of Planck at 
that frequency according to [29j . 

As shown in Fig. 6, while the amplitude of our analyti- 
cal model matches with the residual clustering spectrum 
of point sources in the Planck simulation at multipolcs 
of 10 3 , our analytical model has residual point sources 
that are more clustered than the simulated sources. We 
believe this is due to the finite size of boxes used to sim- 
ulate the source distribution by the Planck team. While 
we fix our point sources to the analytical model as shown 
in Fig. 6, we again capture the uncertainty in the ampli- 
tude with a parameter Ppianck (in this case clustering and 
shot-noise combined as in the case of ACBAR data) and 
marginalize over this parameter when estimating cosmo- 
logical parameters. 

We follow the same procedure as fitting existing 
WMAP and ACBAR data to extract cosmological pa- 
rameters with Planck. We use data at I < 2000, though 
Planck analysis can be extended to higher multipolcs, 
which are likely to be contaminated by additional sec- 
ondary anisotropics beyond SZ_[3l|, and complications of 
the non-Gaussian covariance [32fl . The best-fit cosmolog- 
ical parameters and 68% confidence errors for the stan- 



1 http:/ /www. planck.fr/articlc334. html 

2 http://galics.iap.fr/ 
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FIG. 5: Two-dimensional marginalized distributions showing the 68% and 95% confidence level contours for n s vs. dn s /d\nk 
(left) and n a vs. r (right) of interest with WMAP 5-year and ACBAR data (empty contours) and and WMAP 5-years 
data+ACBAR data with clustering for point sources (blue contours). 



TABLE III: Mean values and marginalized 68% c.l. limits for several cosmological parameters from Planck mock data with 
and without a point source contribution. 



Parameter 


Planck mock data 


Planck mock data 
with clustering 


Planck mock data 
with clustering 




Planck mock data 
with clustering 






Point sources ignored 


Cf x (0 < Ppi anck < 1) 


Shot- 


noise only (Cf™ = Ppianck) 


fl b h 2 


n ri997n+ aoo(:l25 

U.UZZ(U_ 000 24 


0.02753 ± 0.00028 


0.02271 ± 0.00025 




271 + - 00025 

U.UUZ(1_ .00024 


fi c h 2 


0.1082 ±0.0019 


0.0924 ±0.0017 


n in8n+ 0019 
0.1080_ 0020 




in92 +0 0020 
u.iuyz_ 0019 




0.750 ±0.010 


0.831 ±0.007 


n 7 ei +0.011 

U. "1 — 0.010 




0.745 ± 0.010 


n s 


u.yoz_ 007 


1.159 ±0.007 


Qfi1+ 008 

U.yDl_ .007 




0.973 ± 0.007 


T 


0.090 ± 0.007 


260+ 015 


o nqo+ 007 

u.uyu_ .oo6 




n94 +0 007 
u.uy^_ oo8 


A 2 


(2.41 ±0.10) ■ 10" 9 


(1.85 ±0.10) ■ 10 -9 


(2.41 ±0.10) • 10" 9 




(2.40 ±0.10) • 10" 9 




n 7SQ+0-008 
U. r»9_ .009 


0.903±°;SS 


0.787 ±0.009 




0.808 ± 0.009 




0.250 ±0.010 


0.169 ±0.007 


0.2491°;™ 




0.255 ±0.010 


H 


72 4+ 1 - 


84.3±i2 


72.5 ± 1.0 




72.0 ± 1.0 


Zreion 


11.3 ±0.6 


19.9+°-? 


11.3 ±0.6 




11.6 ±0.6 


to 


13.69 ±0.04 


13.04i° ;° 4 


13.69±g : g* 




13.70 ±0.04 


A sz 


1.40 ±0.07 


< 2.00 


imt°oil 




1 84 +aie 

L -°^-0.14 


-Ppianck 






< 1.00(2cr) 


cr 


< 3.6 x 10" 5 ^k 2 -sr (2a) 



dard 6-paramcter ACDM case complemented by Asz and 
-Ppianck are tabulated in Table III. Without point sources, 
we recover the best-fit cosmology that was used to create 
the mock. Once the mock includes point sources and we 
ignore the effect of point sources when model fitting the 
data, we find that the parameters are significantly biased; 
in some parameters this bias is more than 20a. As in the 
case of current data, we consider the possibility that if it 
is adequate to model point sources with just a shot-noise 
power spectrum. We allow Cf n = -Ppianck x 0.0075 /ik 2 -sr 
and fit the data by varying Ppianck- The values of cosmo- 
logical parameters with the shot-noise marginalized over 
is tabulated in Table III, in addition to the 2a upper 



limit on Cf™ from the data. With a shot-noise descrip- 
tion only in the model fit, we find biases in cosmological 
parameters at the level of 1.5er, for example, in the case 
of the spectral index. 

Once we include a model for point source clustering, 
in addition to the shot-noise, and marginalize over the 
overall amplitude with our parameter Ppianck to capture 
the overall clustering amplitude, we find that the biases 
in best-fit parameters from the values used for the mock 
are removed. We show couple of examples for combina- 
tions involving the scalar spectral index n s in Fig. 7 and 
Fig. 8 for Q, c h and as, respectively. As shown, once 
point source clustering is included in the fit, cosmologi- 
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FIG. 6: Planck mock data with and without a contribution 
from point sources: an offset in the temperature power spec- 
trum is visible for small scales. Also showed are the model 
for point sources used in the paper and a simulated signal for 
point sources from public point source maps for Planck HFI 
channels (in this case at 143 GHz) from the Planck team (see 
text for details). 

cal parameter biases are negligible. Ignoring point source 
clustering, however, impacts the measurements signif- 
icantly, though including only a shot-noise for Planck 
point sources still results in an appreciable bias. To 
be completely safe, we suggest that a reasonable model 
for point source clustering and shot-noise be included in 
the cosmological parameter analysis with Planck (at each 
channel used for cosmological measurements) and the un- 
certainty in the modeling or predicting the total point 
source contribution be marginalized over. 

IV. SUMMARY 

The faint radio point sources that are unresolved in 
cosmic microwave background (CMB) anisotropy maps 
arc likely to be a biased tracer of the large-scale structure. 
While the shot-noise contribution to the angular power 
spectrum of radio point sources has been considered so 
far when extracting cosmological parameters with CMB 
data, we have shown here that one should also allow for 
the possibility of source clustering. This is especially nec- 
essary at high frequencies where the clustering of sources 



is expected to dominate the shot-noise level of the angu- 
lar power spectrum at tens of arcminute angular scales. 
As we find, the differences seen by the WMAP team for V 
and W-band angular power spectra do allow point source 
clustering, though one can wrongly conclude clustering is 
unnecessary if lower frequency data are included. 

Here, we have made an estimate of the clustering of un- 
resolved radio sources in both WMAP and ACBAR by 
making use of existing counts at 95 GHz and by making 
several assumptions on the sources such as the redshift 
distribution. To account for the uncertainty in modeling 
the clustering, we included an extra nuisance parameter 
for each dataset and have marginalized over this param- 
eter when model fitting for cosmological parameters. For 
the combination of WMAP 5-year data and ACBAR, we 
find that the spectral index changes from a mean value of 
0.963 ± 0.014 without point-source clustering to a value 
of 0.959 ± 0.014 when the clustering of point sources are 
included in model fits, a difference of 0.3er. We also 
discussed the full parameter set with clustering of radio 
point sources and changes to additional parameters such 
as dn s /dlnk and the tensor-to-scalar ratio r. While we 
find that the differences arc marginal with and without 
source clustering in current data, we have suggested that 
it is necessary to account for source clustering with future 
datasets such as Planck, especially to properly model fit 
anisotropics at arcminute angular scales and using high- 
frequency data. For Planck, we find that simply includ- 
ing the point sources as a shot-noise only out to I of 2000 
for cosmological parameter estimation results in biases at 
the level of 1.5er. While we simply model Planck point 
sources with a single power spectrum, since at high fre- 
quencies both radio and far-IR sources are expected to 
contribute, it may be necessary to return to a proper 
model of total unresolved source clustering in Planck in 
future. 
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